Vacancy ordering and phonon spectrum in the Fe superconductor Ko.sFei 6Se2 
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We have performed Raman-scattering measurements on a high-quality single crystal of the recently 
discovered Fe-based superconductor K0.8Fe1.6Se2 (Tc = 32 K). At least thirteen phonon modes 
were observed in the wave number range 10—300 cm~^. The spectra possess a four-fold symmetry 
indicative of bulk vacancy order in the Fe-deficient planes. We perform a vibration analysis based 
on first-principles calculations, which both confirms the ordered structure and allows a complete 



mode assignment. We observe an anomaly at Tc in the 180 cm 
specific type of electron-phonon coupling. 

PACS numbers: 74.70.-b, 74.25.Kc, 63.20.kd, 78.30.-j 



An mode, which indicates a rather 



I. INTRODUCTION 

Since the discovery of high- Tc Fe-based superconduc- 
tors almost three years ago, much effort has been devoted 
both to exploring their physical properties and to the 
search for new materials in this class. Until recently, five 
such systems had been synthesized and studied, among 
which Fe(Se,Te) is of particular interest as the first Fe- 
based superconductor not requiring the poisonous ele- 
ment As. Although the superconducting transition tem- 
perature (Tc) is low, it has been found for FeSe that a 
maximum Tc of 37 K is obtained under a pressure of 6 
GPa.^ This strong dependence has inspired efforts to find 
new FeSe-based superconductors by introducing internal 
chemical pressure.^ This idea was realized very recently 
in Ko.8Fe2Se2,^ where Tc ~ 31 K is close to that of the 
"122" materials. The system was also identified as being 
isostructural with BaFe2As2^ Several groups have now 
reported similar Tc values by substitution of other alkali 
metals, including Rb and Cs."^ 

The electronic and magnetic properties of the 
Aa;Fe2-2ySe2 compounds have been probed extensively. 
Angle-resolved photoemission spectroscopy (ARPES) 
measurements^^ ^ agree on the presence of one Fermi sur- 
face around the M point, and on a conventional s-wave 
pairing symmetry, but differ on further details. Nu- 
clear magnetic resonance (NMR) measurements in su- 
perconducting crystals find very narrow line widths^ 
strong coupling, singlet superconductivity but with no 
coherence peak, and very weak spin fluctuations. In- 
frared optical conductivity measurements have identified 
the parent compound as a small-gap semiconductor and 
not a Mott insulator. Neutron diffraction studies find 
a very large ordered moment and high magnetic transi- 
tion temperature.^^ High-pressure experiments by several 
groups have led to contradictions which indicate extreme 
sensitivity to the nature of the defects. 

These high-pressure measurements are not the only in- 
dications for a crucial role of Fe vacancies. Superconduc- 
tivity itself occurs only in Fe-deficient samples. A broad 
resistivity peak around 150 K is reduced dramatically by 
altering the Fe content The unexpected abundance of 
infrared- active phonon modes^^ may also be attributed 



to Fe vacancies. Resolving the distribution of Fe vacan- 
cies on the microscopic level is clearly essential for a full 
understanding of the electronic and magnetic properties. 
A study exploiting the unique sensitivity of Raman scat- 
tering to probe local lattice symmetries is thus required. 

Here we measure Raman-scattering spectra in a high- 
quality superconducting crystal of K0.8Fe1.6Se2, as de- 
scribed in Sec. II. Our results, presented in Sec. HI, con- 
tain at least thirteen phonon modes, far more than would 
be expected for a normal 122 structure. The spectra we 
obtain by rotating the crystal show a four-fold symme- 
try, which we identify as C^h oi" C4. These results indi- 
cate an ordered structure of Fe vacancies, and in Sec. IV 
we assume the highest-symmetry defect configuration to 
perform a first-principles vibration analysis. The calcu- 
lated mode frequencies and symmetry assignments are 
fully consistent with our data, confirming completely the 
vacancy ordering pattern. Our temperature-dependent 
spectra, shown in Sec. V, reveal one specific Ag mode 
with a clear frequency jump at Tc, which suggests a close 
connection between phonons and superconductivity. Sec- 
tion VI contains a brief summary. 



II. MATERIAL AND METHOD 

The crystal was grown by the Bridgeman method. 
The detailed growth procedure may be found elsewhere.^- 
The most accurate available composition determina- 
tion is the neutron diffraction refinement, which gives 
K0.8Fe1.6Se2. X-ray diffraction patterns indicate no dis- 
cernible secondary phase. The resistivity was measured 
by a Quantum Design PPMS and the magnetization by 
the PPMS vibrating sample magnetometer (VSM): both 
quantities, shown in the inset of Fig. 1, exhibit transi- 
tions at Tc ~ 32 K, the diamagnetic transition appearing 
particularly sharp. Although the precise demagnetiza- 
tion factor is not available, the estimated superconduct- 
ing volume fraction is (from the sharp diamagnetism and 
absence of a secondary phase) close to 100%. These re- 
sults indicate that the crystal is of very high quality. 

The crystal was cleaved in a glove box and a flat, 
shiny surface was obtained. The surface is parallel to 
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FIG. 1: (Color online) Raman spectra of a K0.8Fe1.6Se2 crys- 
tal at 9 K. Labels d and Cs denote the incident and scat- 
tered light polarizations. Black and red arrows indicate re- 
spectively the assigned Ag and Bg modes, the blue arrow a 
weak mode whose symmetry remains unknown. The accumu- 
lation time for the cross-polarization measurements is approx- 
imately three times that in parallel polarization. The spectra 
labeled Ag and Bg were measured by controlling the relative 
angle between the incident and scattered light polarizations 
once the mode symmetries had been assigned; this assignment 
was performed by combining rotation measurements (Fig. 2) 
with the formal symmetry analysis [Eq. (1)]. Inset: resistivity 
and magnetization data for the same crystal. 



the (a6)-plane, and forms the basal plane for our Raman 
measurements. The two pieces of freshly-cleaved crystal 
were sealed under an argon atmosphere: one was trans- 
ferred into the cryostat within 30 seconds for Raman- 
scattering studies and the other used for symmetry inves- 
tigations by rotating the crystal. The cryostat was evac- 
uated immediately to a work vacuum of approximately 
10~^ mbar. Raman measurements were performed in a 
pseudo-backscattering configuration with a triple-grating 
monochromator (Jobin Yvon T64000), delivering a spec- 
tral resolution better than 0.6 cm~^. The light source 
is a 532 nm solid-state laser (Torus 532, Laser Quan- 
tum) whose beam was focused into a spot of diameter 
ca. 10-20 /im on the sample surface. The beam power is 
lower than 1 mW, and the real temperature at the spot 
was calibrated using the intensity relation between the 
Stokes and anti-Stokes spectra. 



III. PHONON RAMAN SPECTRA: SYMMETRY 
AND VACANCY ORDER 

The phonon modes of K0.8Fe1.6Se2 are shown in the 
low-temperature spectra of Fig. 1. The most strik- 
ing feature is that at least thirteen modes are present 
in the spectra collected with both parallel and cross- 
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FIG. 2: (Color online) Raman spectra obtained by rotation, 
collected in air at room temperature immediately after cleav- 
age. The zero-degree point is chosen arbitrarily (not mea- 
sured relative to the crystal principal axes). Dashed lines 
mark two modes, labeled as ^Bg and ^Ag^ whose relative in- 
tegrated intensity variation is shown as a function of angle 
in the left inset [the dashed line has the form sin^ 20] . Right 
inset: schematic representation of ordered vacancy configura- 
tion with 20% Fe deficiency and four-fold symmetry. 

polarization configurations. All of these modes are 
located below 300 cm~^ (Fig. 4 illustrates this most 
clearly). Our results are therefore consistent both by 
number and location with the infrared optical spectraJ^ 
We stress that only four Raman-active modes are present 
in 122,^^ to which K0.8Fe1.6Se2 was thought to be 
isostructural, and that the spectra of Fig. 1 are com- 
pletely different from those observed in Fe(Se,Te).^^ The 
surprisingly large number of Raman and infrared phonon 
modes reinforces the assumption that the Fe deficiency 
drives a local breaking of lattice symmetry. A vibration 
analysis of the modes in Fig. 1, which we perform below, 
is thus significantly more complicated than for 122 and 
requires additional structural information. 

In Fig. 2 we present spectra obtained by rotating the 
crystal. These display clearly a four-fold symmetry, 
which is highlighted in the right inset. A formal anal- 
ysis based on four- fold symmetry proceeds as follows. 
There are two possible Raman tensors, one for C^hjC^ 
and the other for D^^jD^?^ If d is defined as the angle 
between the crystallographic a-axis and the polarization 
of the incident light, and as the angle between the po- 
larizations of the incident and scattered light, then the 
angle-dependences of the Raman intensities for C^hlC^ 
take the form 

Ia^ = |acos((/p) +csin((/p)p, 

Ib^ = |^^cos((/:) + 2(9) + esin((/9 + 2l9)|^ (1) 

where a, c, and e are matrix elements of the Raman 
tensors, and depend on linear combinations of polariz- 
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FIG. 3: (Color online) Displacement patterns for selected 
Raman-active Ag [(a) and (b), labeled as ^Ag and ^Ag in Ta- 
ble I] and Bg modes [(c) and (d), labeled as "^Bg and ^Bg in 
Table I] of K0.8Fe1.6Se2. Fe atoms connected by red lines have 
right-handed chirality. Amplitude arrows in (d) are scaled up 
by a factor of 1.4 over the other panels. Atomic structures 
and displacement patterns prepared using XCRYSDEN.— 



ability derivatives for the atoms involved in each mode. 
Similarly, for D^h/D^ one has 



Ia,, = |acos((/:))p, 



|esin((/:) + 2(9)p. (2) 



Because rotating the crystal varies the angle 0^ it is clear 
that modes with constant intensities are all A- type {Aig^ 
Ag, Ai, or A), while those showing an intensity mod- 
ulation are 5-type. Further, Eq. (2) prohibits the ob- 
servation of A-type modes in D^h symmetry in a cross- 
polarized configuration {(p = 7r/2), which contradicts the 
observations in Fig. 1. Thus our results exclude explicitly 
the D4/J, symmetry expected for the 122 structure, and 
allow us to deduce that the crystal symmetry is C4/1 or 
C4. The A- type modes in Fig. 1 may therefore be identi- 
fied as Ag or A phonons. This symmetry constrains very 
strongly the possible vacancy ordering patterns. 

Indications for vacancy ordering in Aa;Fe2-2ySe2 sys- 
tems have been obtained in several recent measurements. 
Imaging by transmission electron microscopy (TEM)^^ 
reveals rather well-defined surface order for y = 0.5 and 
y = 0.2. NMR measurements^^ find additional line- 
splittings. While Raman scattering is a surface tech- 
nique, the wavelength, laser spot size (10-20 /im), and 
penetration depth mean that macroscopic symmetries are 
probed. The 20% vacancy concentration in our crystal 
{y = 0.4) suggests the structure illustrated in the inset 
of Fig. 2: this arrangement has a ^/bx^/b unit cell, has 



local four-fold symmetry, and, crucially, has lost the two 
mirror planes perpendicular to the (a6)-plane, in com- 
plete consistency with our Raman measurements. This 
microscopic structure has been demonstrated spectacu- 
larly bv the very recent neutron diffraction studies of 
Ref. [l2|- X-ray diffraction studies of several crystals with 
y = 0.37-0.39 have confirmed this result.— 



IV. FIRST-PRINCIPLES CALCULATIONS 

To compare with the observed phonon modes, we 
have calculated the nonmagnetic electronic structure and 
zone-center phonons of K0.8Fe1.6Se2 from first-principles 
density functional calculations. We use the Vienna ab- 
initio simulation package, which is based on the projec- 
tor augmented wave (PAW) method^ ^ with a general gra- 
dient approximation (GGA)^^ for exchange-correlation 
potentials. The structural information determined by 
neutron diffraction^ is that K0.8Fe1.6Se2 has space group 
I4/m and point group C^h- The Wyckoff positions of the 
K, Fe, and Se atoms and the corresponding symmetry 
analysis are listed in Table I. The frequencies and dis- 
placement patterns of the phonon modes were calculated 
using the dynamical matrix method. 

In total, 17 Ag and Bg modes are allowed by the crys- 
tal symmetry, among which 11 of the modes in the po- 
larized Raman spectra of Fig. 1 may be assigned accu- 
rately from the calculations (Table I). The vibration pat- 
terns of selected modes are shown in Fig. 3. The cal- 



TABLE I: Symmetry analysis for I4/m and assignment of 
Raman- active vibration modes in K0.8Fe1.6Se2 
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FIG. 4: (Color online) (a) Raman spectra at selected tem- 
peratures. An example of fitting with Lorentz functions is 
shown at the base of the panel for 300 K data collected in 
parallel polarization with an accumulation time four times 
that of the other spectra. The frequency resolution is then 
0.3 cm~^, which allows the explicit separation of 13 modes in 
the narrow frequency interval 50-300 cm~^. The fits used to 
obtain the data points in Figs. 4(b) and 5 have a coefficient 
of determination > 0.95. (b) Temperature-dependence of 
the relative intensities of the 66 cm~^ and 274 cm~^ modes. 
The solid line shows the Bose-Einstein thermal factor. 



culated phonon frequencies, 75.1 and 212.6 cm ^ for Ag 
modes [Fig. 3(a) and (b)] and 238.3 and 279.0 cm-^ for 
Bg modes [Fig. 3(c) and (d)], agree very well with the 
66.3, 202.9, 214.3, and 274.9 cm"^ modes in the experi- 
mental measurements. The Ag and Bg modes shown in 
Figs. 3(b) and (c) correspond exactly to the Aig mode 
of As and Big mode of Fe in the 122 system. The suc- 
cess of the phonon assignment confirms completely the 
assumed pattern of vacancy ordering. Only three of the 
observed 13 modes cannot be assigned well, even though 
two have Ag character. We suggest that these may origi- 
nate from distortions of the Fe-Se layer, which lie outside 
the symmetry constraints applied. 



V. TEMPERATURE-DEPENDENT SPECTRA 

We have also measured Raman spectra over the 
full temperature range to 300 K [Fig. 4(a)]. The 
temperature-dependence of the phonon modes is de- 
scribed explicitly by fitting with Lorentz functions. All 
relative intensities show a conventional phonon form, 
with one exception (below). This line-shape analysis 
shows no evidence of the Fano asymmetry character- 
istic of electron-phonon coupling. However, all of the 
phonons are rather broad when compared with those in 
122 systems, due perhaps to intrinsic disorder (incom- 
pletely ordered vacancies). This analysis alone cannot 
therefore exclude a generic electron-phonon coupling in 
K0.8Fe1.6Se2. 

The Ag mode at 66 cm~^ behaves rather anomalously. 
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FIG. 5: (Color online) (a)-(d) Temperature-dependence of 
phonon frequencies for four modes with the highest measured 
intensities. Insets magnify the region around Tc. Horizontal 
and vertical scales are the same in all panels and insets. 



The steep fall in its relative intensity [Fig. 4(b)] fol- 
lows a Bose-Einstein form only to lowest order, but with 
significant deviations. In addition to a clear intensity 
jump at Tc, there is a further apparent anomaly around 
160 K, where X-ray diffraction studies find no structural 
transitions^ We have confirmed (Table I) that the 66 
cm~^ phonon is the vibration mode of a Se atom, and 
suggest that this is the origin of its unusual behavior. 

The temperature-dependence of selected phonon fre- 
quencies is shown in Fig. 5. The superconducting transi- 
tion has very little effect on the frequency in every case 
but one: the Ag mode at 180 cm~^ exhibits a clear jump 
of approximately 1 cm~^ at Tc [Fig. 5(a)]. This anomaly 
is definite evidence for a particular type of connection be- 
tween phonons and superconductivity, and it constitutes 
a statement concerning both symmetry and energies. For 
the former, the symmetry relation between the super- 
conducting gaps on the available Fermi surfaces should 
be consistent with Ag. For the latter, we suggest that 
a couphng of this phonon mode to the superconducting 
quasiparticles may be visible^^ in ARPES spectra, where 
it should be expected at energies E ~20— 25 meV and at 
the r point (the part of the Brillouin zone probed by Ra- 
man scattering). Anomalous features are indeed present 
in this region in the two most recent ARPES studies of 
Aa^Fe2-2ySe2 superconductors,^ and we propose that high- 
resolution ARPES measurements may reveal the physics 
behind this phenomenon. 



VI. SUMMARY 

To conclude, we have performed Raman-scattering 
measurements on the recently discovered Fe-based su- 
perconductor K0.8Fe1.6Se2. Using a high-quality single 
crystal, we find that there exist at least thirteen phonon 
modes and that the spectra have four-fold symmetry. 
This demonstrates the presence of long-range-ordered 
configurations of Fe vacancies, which dictate the local lat- 
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tice vibrations. We perform first-principles calculations 
to obtain a complete and consistent phonon assignment, 
confirming the nature of the vacancy ordering pattern. 
We find that only one mode exhibits a change in fre- 
quency around Tc, suggesting a rather specific connection 
between superconductivity and lattice vibrations. 
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